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The maximum-likelihood principle unifies inference of quantum states and processes from exper- 
imental noisy data. Particularly, a generic quantum process may be estimated simultaneously with 
unknown quantum probe states provided that measurements on probe and transformed probe states 
are available. Drawbacks of various approximate treatments are considered. 
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I. INTRODUCTION 

The various quantum-state reconstruction techniques 
developed during recent years have made it possible to 
completely reconstruct an unknown state of a quantum 
mechanical system provided that many identical copies 
of the state are available. These reconstruction methods 
are nowadays routinely applied to the evaluation of the 
experiments where quantum states are generated, ma- 
nipulated and transmitted. The field was pioneered in 
the beginning of nineties in quantum optics, where the 
optical homodyne tomography has been devised for re- 
construction of the quantum state of traveling light field 
El |2[ ||, ||, |). Since then, many other reconstruction 
methods applicable to various physical systems have been 
developed §, |, |, | 0, £§ || g g || 
The inference of quantum states plays very important 
role in the present-day experiments |p], |2(], [2l], |2^] . 

Most of the reconstruction methods, such as the direct 
sampling in optical homodyne tomography, are based on 
a direct linear inversion of the experimental data. This 
approach is conceptually simple and feasible. However, it 
may lead to certain unphysical artifacts such as the neg- 
ative eigenvalues of the reconstructed density matrix. In 
order to avoid these unphysical artifacts, an estimation 
method based on statistical maximum-likelihood princi- 
ple has been devised for the reconstruction of a generic 

quantum state || HI fH III & This a PP roach 
guarantees the positive semidefiniteness and trace nor- 
malization of the reconstructed density matrix. These 
necessary conditions are incorporated as constraints, so 
as a certain prior information from the statistical point 
of view. Remarkably, the maximum likelihood estima- 
tion can be interpreted as a genuine generalized quantum 
and can be related to the infor- 
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measurement 

mation gained by optimal measurement and the Fisher 
information p9| . 

Given current interest in the quantum-information 
processing, it is of paramount importance to reconstruct 
not only the quantum states but also the transforma- 
tions of these states — the quantum mechanical processes. 
The examination of quantum communication channels 
and the evaluation of the performance of quantum gates 
are the examples of practical applicability of quantum- 
process reconstruction [|o[ ||l], |3^, [33], [34|, |35|, |36|, [37| . 



All necessary properties of the deterministic quantum 
transformations, namely the complete positivity and 
trace preservation can be again incorporated within the 
maximum-likelihood approach as the appropriate con- 
straints |38| ] . Compared with other reconstruction meth- 
ods the maximum-likelihood approach seems to be com- 
putationally more difficult. Therefore several simplifi- 
cations and approximations of the maximum-likelihood 
technique have been suggested recently |39[ [|0| . 

In this paper we present a unified approach to the 
maximum-likelihood reconstruction of quantum states 
and quantum processes. Extremal equations for the re- 
constructed quantum state and for quantum process are 
derived in Section II. These equations can easily be solved 
numerically by means of repeated iterations. Particular 
attention will be paid to the probing of the quantum 
process by entangled states which attracted considerable 
attention recently. In Section III we consider a realis- 
tic scenario where an unknown quantum transformation 
is probed by unknown states and the measurements are 
performed on both the input and output states. We pro- 
pose a method for simultaneous estimation of the un- 
known probe states and the quantum process from the 
collected experimental data. The comparison of the exact 
maximum-likelihood method with the approximate ones 
is carried out in Section Finally, the conclusions are 
given in Section V. 



II. RECONSTRUCTION OF QUANTUM 
PROCESS 



Let us start with a brief review of the maximum- 
likelihood reconstruction of a quantum state. We assume 
a finite number N of identical samples of the physical sys- 
tem, each in the same but unknown quantum state de- 
scribed by the density operator p. Having these systems 
our task is to infer the unknown quantum state p from 
the results of the measurements performed on them. We 
consider the positive operator-valued measure (POVM) 
pif that yields probabilities pi of individual outcomes, 



pi = Tr [pIli] , Pi>0, ^2pi = 1 - 



(1) 



If the POVM n; is tomographically complete it is possi- 
ble to determine the true state p directly by inverting the 
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linear relations ([j]) between the probabilities pi and the 
elements of the density matrix p. However, there is no 
way how to find out the exact probabilities p\ since only 
a finite number N of samples of physical systems can be 
investigated. In the case of Ni occurrences of outcomes 
11; the relative detection frequencies fi — Ni/N represent 
the only data that could be used for reconstructing the 
true state p. The maximum-likelihood approach to this 
reconstruction problem consists in finding a density op- 
erator p est that generates through Eq. (|l]) probabilities 
Pi which are as close to the observed frequencies / ; as 
possible ||,|2f], 



Pest = arg max C[fi,pi(p)}, 
p 

£[fl,Pl(p)} = 



(2) 
(3) 



The measure £[fi,Pi(p)] of the distance between the 
probability distribution pi and the detected relative fre- 
quencies // seems to be arbitrary. However, it can be 
shown that the reconstruction procedure can be inter- 
preted as a generalized POVM measurement if the log- 
likclihood measure (§) is used @ JJ6|. The maximum- 
likelihood principle has been successfully applied to many 
problems of quantum-information processing, for exam- 
ple to reconstruction of the spin state of an electron or po- 
larization state of a photon p5[ , reconstruction of entan- 
gled spin state J28| ] , estimation of quantum measurement 
|t2j , design of the optimal discrimination device for com- 
munication through a noisy quantum channel Eland 
characterization of the universal cloning machine [B4J . 

The challenging problem of the maximization (0) of 
the log-likelihood functional (||) on the space of positive 
semidefinite operators p, Tr[p] = 1, has been treated with 
the help of the numerical up-hill simplex method |p7f . A 
more analytical approach to the problem involves a for- 
mulation of nonlinear extremal operator equation for the 
density matrix that maximizes the log-likelihood func- 
tional [E3l uA ESl , 



P = p- 1 Rp, R = J2- Il l, 

where the Lagrange multiplier p, reads 
M = Tr[ J Rp]=^/ i = l. 



(4) 



(5) 



The crucial advantage of the equation (Q) is that it is 
suitable for iterative solution, as has been demonstrated 
on many particular reconstruction problems. A combi- 
nation of equation (0) and hermitian conjugate equation 
leads to the symmetric extremal equations in the mani- 
festly positive semidefinite form |}42|, 



(6) 



p = p~ 2 RpR, p = (Ti-[RpR}) 1/2 . 
The iterations 

p {n+l) =IJ {n)- 2 R {n) p {n) R {n) (J) 



preserve the positive semidefiniteness and trace normal- 
ization of the density operator p. 

While density operator describes the state of phys- 
ical system, the linear completely positive (CP) map 
describes the generic transformation of physical system 
from quantum state p m to quantum state p ou t- The 
mathematical formulation of CP maps relies on the iso- 
morphism between linear CP maps Ais from operators 
on the Hilbert space TC to operators on the Hilbert space 
K, and positive semidefinite operators S on Hilbert space 
H®K |45l pi Wft. 



Pont 



= Ms[pin]=Tr n [Sfyl K ] 



(8) 



where TLjc is an identity operator on the space K. and 
T denotes the transposition. The deterministic quantum 
transformations preserve the trace of the transformed op- 
erators, Tr^[/9 out ] = Tr«[pi n ]. Since this must hold for 
any p- m the operator S must satisfy the condition 



Tr K [S\ = l n , 



(9) 



where TL-h is an identity operator on space H. The condi- 
tion (g) effectively represents (dimW) 2 real constraints. 

Making use of the formalism (g) we may formulate 
the exact maximum-likelihood principle for estimated CP 
map S in a particularly simple and transparent form and 
we can also straightforwardly extend the results obtained 
in Ref. |38| to the cases when the input and output 
Hilbert spaces have different dimensions. 

Let p m denote the various input states from the space 
Ti. that are used for the determination of the quantum 
process. Measurements described by POVMs U m i are 
carried out on each corresponding output state from 
space /C. Let / m ; denote the relative frequency of de- 
tection of the POVM element H m i- The estimated op- 
erator S should maximize the constrained log-likelihood 
functional || || 

£c[fmUPml{S)] = ^ f ml ^ n Pml - Tf[A5], (10) 



1 n . I 



Pn 



Tr [S pi 



(11) 



where A = A ® Ij; and A is the matrix of Lagrange mul- 
tipliers that account for the trace-preservation condition 
(g). The extremal equations for S can be obtained by 
varying functional (fiXm with respect to S, which leads to 



S = A KS, 



r , \ ^ fml T 

, Pral 



n 



ml ■ 



(12) 



Further we have from Eq. (|T^) that S = SKA^ 1 . When 
we insert this expression in the right-hand side of Eq. 
( |l2| ) , we finally arrive at symmetrical expression suitable 
for iterations, 



S = A KSKA~ 



(13) 



The Lagrange multiplier A must be determined from the 
constraint (g). On tracing Eq. @ over space K, we 
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obtain quadratic equation for A which may be solved as 
\ = (Tr K [KSK]) 1/2 . (14) 

The operator A is positive definite because KSK is 
positive definite operator. The system of coupled Eqs. 
( |l3| ) and ([l4]) may be conveniently solved numerically by 
means of repeated iterations, starting from some unbi- 
ased CP map, for example = l^®;c/(dim/C). It is 
important to note that Eq. ( |l3| ) preserves the positive 
scmidefinitcness of S and also the constraint Tr/cfS] = 
t-H is satisfied at each iteration step. 

The density matrix S representing the CP map Ms 
can be in fact prepared physically in the laboratory if we 
first prepare a maximally entangled state on the Hilbcrt 
space Ti(B)Ti and then apply a CP map to one part of 
this entangled state. In this way the quantum-process 
tomography can be transformed to the quantum-state 
tomography. More generally, this suggests that it may 
be useful to employ entangled quantum states as probes 
of the unknown quantum process [ |37| . 

Let p m ,AB denote the entangled state on the Hilbert 
space Ha ®Wb that serves as a probe of the CP map S 
that is applied to the subsystem A. A joint generalized 
measurement described by the POVMs H m i if performed 
on the output Hilbert space K ® Hb ■ The log-likelihood 
functional has the form (|l0|), only the formula for the 
probability p m \ changes to 



Pml 



= Tr 



HaHbK 



[(S<8l Ws )(p 



T A 



,AB 



1jc)(1 



IT, 



(15) 

where Ta stands for the partial transposition in the sub- 
system A. Consequently, the operator K appearing in 
the extremal Eqs. ( |l3| ) and (fli) must be calculated as 
follows, 



K 



Y,—Trn B l( P l A AB ®^Kln A ®Il m i)}- (16) 

, Pml 



Apart from these modifications of p m ; and K one can 
proceed as before and solve Eqs. (jl^) and (|lj) by means 
of repeated iterations. 



III. QUANTUM PROCESS MEASUREMENT 
BY UNKNOWN PROBE QUANTUM STATES 

Up to now quantum states and processes have been 
treated independently. However, this is just a simpli- 
fication typical for the realm of physical experiments. 
Widely accepted strategy how to approach a complex 
problem is to specify some partial subproblems, address 
them separately and merge the solutions. This technique 
usually gives good answer in the technical sense. Though 
this is possible even in quantum theory, there are no fun- 
damental reasons for such a factorization. To consider 
the full problem without splitting it into isolated sub- 
problems is technically more advanced but could be ad- 
vantageous. This strategy will be demonstrated on the 
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FIG. 1: Scheme of setup for the generalized measurement of 
quantum process using unknown quantum states as probes. 



synthesis of the problems treated separately in the previ- 
ous section. Let us assume the estimation of the generic 
process with the help of set of probe states, identity of 
which is also unknown. What is only known to the ex- 
perimentalists are the output of certain measurements 
performed on the ensemble of probe states and on the 
ensemble of transformed probe states. In this sense all 
the considerations are done ab initio, since only results 
of generic measurements are required. A quantum object 
could be considered as known only to the extent speci- 
fied by some preceding measurements. All the physically 
relevant results will be derived exclusively from the ac- 
quired data, where input states and their transformation 
are inseparably involved. States and their transforma- 
tion should be considered as quantum objects. As such 
they are affected by quantum fluctuations, since in every 
experiment a certain portion of the noise will be present 
on the microscopic level. 

In the following the probe quantum states p m will be 
treated as unknown mixed states and they will be inferred 
together with the unknown quantum process S. In accor- 
dance with the theory presented above let us consider the 
set of probe states p m on the space H. By means of un- 
known quantum process S these states are transformed 
onto output states p m ,out in the space fC. The observa- 
tion must be more complex now involving the detection 
on the ensemble of both the input and the output states. 
For this purpose the corresponding POVM elements will 
be denoted by ir m k and n m ;. The diagram involving de- 
tected signals and measurements is shown in Fig. ^. Let 
f m k denotes the relative frequency of detection of the 
POVM element ir m k in the input space H and F m i de- 
notes the relative frequency of detection of the POVM el- 
ement n m ; in the output space K. The frequencies f m k, 



Sfe fn 



1, and F m i, J2i Fmi = L approximate the 



true probabilities p m k and P m i of individual outcomes, 
respectively, 



Pmk = Tr w [p m 7r m/ ] , 
Pml = Trx;[p m ,outn m; ] = Tr[S"(p^ ® U m i)] 



(17) 



where the relation (||) was used. The estimated process 
S and probe states p m should maximize the constrained 
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log-likelihood functional 



£c = fmk hi jW + F m i In P„ 



(18) 



The additivity of log likelihood reflects the independence 
of observations performed on the input and output states 
with the same degree of credibility. The Lagrange multi- 
pliers [i m and A = A ® lx fix necessary constraints — the 
trace normalization of the states, Tr[p m ] = 1, and the 
trace-preserving property (JsJ) of the process S. 

The coupled extremal equations for the probe states 
p m and for the process S can be obtained by varying JIs| ) 
with respect to independent variables p m and S, which 
yields 



A^KSKA' 1 = S, 



(19) 
(20) 



where 



Efmk 
^mk~ 

, Pmk 
k 



+Tr K 



\ , "ml 



in, I 



■F ml t 
~p~Pm 



(21) 



(22) 



and T-yi is operator of partial transposition in space Ti. 
acting on space TL®JC. The Lagrange multipliers can be 
determined from the appropriate constraints, 



Mm — {Tr?{{R m p m R m ]) - 
A = (Tr K [KSK])* . 



(23) 
(24) 



All necessary properties of the quantum states p m and 
the quantum process S are satisfied during the iterative 
solution of the extremal equations (|l9|)-(p4|). 

In the rest of this section we illustrate the developed 
method on the estimation of a quantum process S that 
transforms one qubit state to another one, dim7i = 
dim/C = 2. The process S under consideration consists of 
a unitary 7r/4- rotation in irz-plane of the Bloch space and 
a subsequent non-unitary damping. The unitary part of 
the process can be represented by its action on ortogonal 
states |0) and |1), 



|0) -> cos6»|0) +sin6»|l), 
1) -> cos6»|l) - sine jo), 



(25) 



where 9 = ir/8. The non- unitary part of the process 
is described by the operator dD + (1 — d) E, where we 
chose d = 1/2. The process D = l-^ 8 ^/2 is totally 
depolarizing channel that maps all states to the maxi- 
mally mixed state and E is the identity transformation. 



We have performed numerical simulations of the i: m k and 
W m i measurements for M = 20 input probe states p m and 
the corresponding transformed states p m ,out respectively. 
The mixed states p m have been randomly generated. We 
consider a convenient experimental realization where the 
same measurements are performed on all input as well 
as output states. In the present example this POVM 
measurement consists of tomographically complete set of 
projective measurements in x, y and z directions, each 
made on N = 1000 identical samples of the probe states 
p m before and after transformation. Therefore, the to- 
tal number 6MN of the probe states have been used up. 
Theoretical probabilities p m k and P m i have been evalu- 
ated according to Eq. (p7|). They represent mean values 
of the multinomial distributions of the relative frequen- 
cies f m k and F m i- Corresponding variances are approxi- 
mately given by p mk {l -p m k)/N and P m i(l - P m i)/N, 
respectively. The data f m k and F m i have been obtained 
by means of Monte-Carlo simulation. Subsequently, we 
have iteratively solved extremal equations (|l9])-(|24|). Re- 
sult of this procedure is shown in Fig. ^. Only 12 real 
independent elements of estimated process are plotted in 
form of a vector {S n }^ =1 . The estimated values are well 
corresponding to the true ones. 
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FIG. 2: Elements of the reconstructed quantum process 
(solid) are compared with the theoretical ones (dashed) for 
the rotating-damping channel and 20 various probe states. 

The simultaneous reconstruction discussed above 
yields a higher likelihood of estimated quantum objects 
than separate reconstructions of probe states and a quan- 
tum process. This seems to be a general rule. The likeli- 
hood £ s i m obtained by simultaneous reconstruction jl8|)- 
(^4~| ) of the quantum process S and the probe states p rn is 
always higher than the sum £ scq of likelihoods C Pm ob- 
tained by the separate reconstructions (||) of the probe 
states p m from data f m k , F m i and likelihood £$ of the 
estimated quantum process S (p^)-(p^), where the recon- 
structed probe states are utilized. The ratio £ S i m /£ seq 
averaged over an ensemble of possible experimental data 
is plotted in Fig. |3| for several numbers of probe states and 
various numbers N of measurements. The true process S 
and the POVM measurements 7r m fe, H m i are the same as 
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in the previous example. A significant improvement is ob- 
tained by using the proposed simultaneous reconstruction 
method in the case of small number N of measurements, 
so in the case of noisy data. The quantitative difference 
between simultaneous and sequential reconstruction pro- 
cedures changes to qualitative one for a tomographically 
incomplete POVM measurement in the input or output 
space. Data acquired by such a measurement could be 
insufficient for the sequential reconstructions, however, 
they can be sufficient for the simultaneous one. For ex- 
ample, projective measurements in x, y directions in the 
input space and projective measurements in y, z direc- 
tions in the output space represent this case. Thus the 
presented simultaneous reconstruction technique is appli- 
cable to the problems, where routine sequential methods 
fail. 

1.12 
1.1 



Accordingly, the exact likelihood functional (||) can be 
replaced by the approximate one J39fl, 




50 100 150 200 250 300 

N 

FIG. 3: The average ratio of the likelihood attained by simul- 
taneous reconstruction of a quantum process and probe states 
and the likelihood attained by sequential one. The process is 
probed by 15 (diamond), 30 (plus) and 45 (square) quantum 
states. The ratio decreases with the increasing number N of 
measurements. 



IV. APPROXIMATE METHODS 

Recently, approximate reconstruction methods based 
on the maximum likelihood have been presented. Two 
ways can be followed to modify the exact maximum- 
likelihood principle — either simplification of the distance 
measure (|3|) [ p9[ or releasing some constraints on quan- 
tum states and processes [|40| . 

For large number N of identical samples of quantum 
states available for inspection going before the state re- 
construction the relative frequencies // fluctuate around 
the true values pi according to the multidimensional 
Gaussian distribution that approximates the exact multi- 
nomial one, 



exp 



(fi-pi) 2 
2a? 



(26) 



(27) 



The reconstruction based on this functional loses the 
essence of the generalized measurement, nevertheless, it 
preserves all physical properties of estimated quantum 
states. The Gaussian limit of the likelihood method 
have been recently applied to the reconstruction of 
polarization-entangled states of light (i^, ^9, M. Un- 
like this, the approximate reconstruction of quantum pro- 
cesses proposed in Ref. [|IJ uses the exact likelihood func- 
tional (|l0|), however, it decreases the number of the con- 
straints incorporated by the Lagrange multipliers. The 
(dirnW) 2 necessary conditions that guarantee the correct 
normalization of the estimated process are replaced by 
a single condition, Tr[S] = dim7Y. This is equivalent to 
assuming that the Lagrange multiplier A is proportional 
to identity operator. 

In order to compare explicitly the exact maximum- 
likelihood estimation of quantum process |3q| with ap- 
proximate method presented in Refs. }44| we have 
carried out extensive numerical simulations. Quantita- 
tive comparison of the two approaches was based on the 
variances of estimates Se (exact) and S*a (approximate), 



4 = <Tr[(S E - 5 trac ) 2 ]) c 
a| = (Tr[(S A -S tl . uc ) 2 ]) c 



(28) 



■pi{l-pi)/N. 



where (. . .) cns denotes averaging over an ensemble of all 
possible experimental data and St rU c denotes the true CP 
map. For a given fixed CP map, input states, and output 
measurements, we have repeated 1000 times a simulation 
of the measurements and reconstruction of the CP maps 
Se and Sa- Subsequently we have calculated variances 
( p8| ) as statistical averages over the acquired ensemble. 
We have found that the exact maximum-likelihood esti- 
mation yields in all cases much lower variance than ap- 
proximate approach. This is a direct consequence of the 
fact that the exact treatment takes into account all con- 
straints imposed by quantum mechanical laws on the es- 
timated operator S. A typical example is shown in Fig. ^|. 
In this case, the quantum process is a unitary transfor- 
mation (|25| ) of a single qubit. Six different input states 
are considered — eigenstates of three Pauli matrices <r x , 
<7y, and a z . 3N copies of each input state are used. On 
each corresponding output state, a spin projection along 
axes x, y and z is measured N times. As can be seen 
in Fig. ^ the variance a E is approximately twice smaller 
than variance a A , which is a significant difference. In 
fact, for CP maps which do not represent unitary trans- 
formations, such as Pauli damping channel, the difference 
may be even stronger. 
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50 100 150 200 250 300 

N 

FIG. 4: The variance erf (diamond) of the exact maximum- 
likelihood reconstruction of a quantum process and the vari- 
ance a\ (plus) of the approximate one for various numbers N 
of measurements. 



CONCLUSIONS 



handle data from many different experimental configu- 
rations such as the probing of quantum processes with 
entangled states or a simultaneous reconstruction of an 
unknown process and unknown states that are used to 
probe this process. The extremal equations (^), (jlq)- 
(|i~4|), and (|l9|)-(|24|) for the most likely quantum state 
and process can be very efficiently solved numerically 
by means of repeated iterations. The exact maximum 
likelihood estimation of quantum objects has been com- 
pared with the approximate methods. The approximate 
ones yield estimates whose variance is typically substan- 
tially larger than in the case of the exact approach. This 
comparison clearly illustrates the importance of keeping 
all the constraints imposed by quantum theory. Loosely 
speaking there is always a choice — either to acquire less 
portion of the data and then to adopt more sophisti- 
cated algorithm for its evaluation or vice versa. The 
efficient and precise reconstruction technique discussed 
in the present paper can find applications in design and 
evaluation of quantum-information devices and contem- 
porary quantum experiments. 



The unified approach to inference of quantum states 
and quantum processes from experimental noisy data 
has been presented. The proposed technique based on 
the maximum-likelihood principle preserves all proper- 
ties of the states and the processes imposed by quan- 
tum mechanics. This method is very versatile and can 
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